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ABSTRACT 


'  ^The  production  of  high  current  density  ion  beams  by  the  direct 
interaction  of  an  intense  electric  field  on  a  liquid  surface  was  investi¬ 
gated.  Electrohydrodynamic  spraying  techniques,  which  have  been  tradi¬ 
tionally  utilized  for  the  production  of  mul t i -molecula r  charged  particle 
beams,  were  optimized  for  the  production  of  ionic  particles. 

Simultaneous  ion  formation  and  acceleration  was  achieved  by  applying 
a  high  DC  voltage  between  a  hemispherical  liquid  meniscus  supported  at  the 
tip  of  a  metallic  capillary  and  a  concentric,  planer  accelerating  electrode. 
With  proper  feed-back  control  of  the  applied  voltage  and  proper  wetting  of 
the  capillary  by  the  working  fluid  the  hemispherical  meniscus  could  be 
deformed  by  the  action  of  tne  applied  electric  field  into  a  stable  conical 
geometry.  At  the  apex  of  this  cone  the  electric  field  intensity  was  of 
sufficient  amplitude  to  field  evaporate  ions  from  the  liquid  surface. 

Such  a  field  evaporation  source  using  a  liquid  anode  can  operate  inde¬ 
finitely  at  high  emitted  current  levels. 

The  field  evaporated  particles  were  analyzed  to  determine  the  useful¬ 
ness  of  this  technique  as  a  source  of  charged  particles  for  electrostatic 
microthrusters.  The  total  emitted  current,  specific  charge  distribution 
of  the  emitted  particle,  and  formation  energy  of  the  particles  were 
measured  in  terms  of  the  operating  parameters  of  the  system. 

Stable  ion  currents  of  up  to  500  microamps  were  obtained  from  a 
single  capillaiy  emitter.  The  ion  beam  consisted  of  75%  gallium  dimers, 

24%  indium  dimers  and  1%  gallium  tetramers.  The  resulting  thrust  produced 
was  73.8  micronewtons. 
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1.  ELECTROSTATIC  PROPULSION 


1 . 1  Mot i va  t ion 

Ever  since  Tsiolkovsk i  i  wrote  his  now  famous  rocket  equation  that 
relates  the  terminal  velocity  of  a  space  vehicle  to  its  exhaust  velocity 
and  mass  ratio,  rocket  designers  have  recognized  the  need  for  a  thrustor 
with  a  high  exhaust  velocity  and  low  propellant  mass  flow  rate.  The 
exhaust  velocity  of  a  chemical  rocket  is  limited  by  the  enthalpy  of  the 
reaction  between  the  fuel  and  oxidi2er.  The  most  energetic  propellant 
combinations  such  as  hydrogen-oxygen  and  hydrogen-fluorine  do  not  yield 
sufficiently  higr.  exhaust  velocities  for  optimum  payload  ratios  or  terminal 
velocities.  Consequently,  a  new  source  of  thrust  was  needed  that  would 
operate  more  efficiently  than  chemical  systems. 

The  use  of  electrostatically  accelerated  ions  to  produce  thrust  was 

(2) 

suggested  by  Goddard  in  1906.  An  electrostatic  thrustor  represents, 
at  least  ideally,  a  highly  efficient  means  of  converting  stored  energy 
into  thrust.  The  thrust  produced  is  usually  very  small  compared  to  that 
of  chemical  rockets,  and  for  this  reason  electrically  propelled  vehicles 
cannot  operate  under  conditions  where  appreciable  atmospheric  drag  is 
encountered  or  where  the  trajectory  demands  propulsion  against  large  gravity 
forces.  Electric  thrustors  function  by  applying  a  small  thrust  for  a  very 
long  period  of  time,  whereas  chemical  rockets  produce  a  very  large  thrust 
for  a  short  time.  For  many  missions  planned  for  the  near  future,  electric 
propulsion  systems  show  marked  advantages  in  overall  operational  efficiency 
over  chemical  rockets^^ 


1 


2 


It  has  been  the  intent  of  this  research  to  further  the  state-of-the- 
art  of  electrostatic  propulsion  systems.  In  the  following  section  the 
physics  of  electrostatic  thrustors  will  be  outlined  and  the  important 
criteria  of  a  beam  of  charged  particles  used  for  thrust  production  will 
be  discussed.  On  the  basis  of  these  criteria,  existing  electrostatic 
propulsion  systems  will  be  explained  nnd  compared.  Finally,  noting  the 
operation  of  current  electrostatic  thirstors  and  the  design  parameters 
imposed  by  the  physical  laws  governing  the  motion  of  space  vehicles,  the 
objectives  of  this  study  will  be  put  forth. 


1 .2  Physics  of  Electrostatic  Thrustors 
A .  Electrostatic  Acceleration 

The  motion  of  a  space  vehicle  is  governed  by  Newton's  laws.  An 
isolated  space  vehicle  can  be  accelerated  by  emitting  a  propellant  of  mass 
m  with  an  exhaust  velocity  v.  The  force  which  expels  the  propellant 
particles  is  equal  to  the  rate  of  change  of  momentum  of  the  propellant 
which  Is  given  by  Newton's  second  law. 


F 


d  (mv) 

dt 


(1-1) 


For  electrostatic  systems  the  exhaust  velocity  v  is  constant  and  is 

determined  by  the  accelerating  voltage  V  and  the  specific  charge  of  the 

a 

propellant  particles  q  through  the  relationship 

s 

v^  =  2V  q  (1-2) 

as 

The  exhaust  velocity  is  plotted  as  a  function  of  the  accelerating 


voltage  for  various  particle  specific  charges  in  Fig.  (1-1). 


FIG.  1-1  EXHAUST  VELOCITY  AS  A  FUNCTION  OF  ACC! 

WITH  PARTICLE  SPECIFIC  CHARGE  AS  A  PAI 
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For  e  constant  exhaust  velocity  system  Newton's  second  law  takes  the 


F  =  mv 


where  m  is  the  time  rate  of  change  of  propellant  mass  or  mass  flow  rate. 
The  same  force  F  that  expels  the  propellant  also  acts  upon  the  vehicle 
as  expressed  by  Newton's  third  law. 

The  thrust  can  be  written  in  terms  of  the  total  emitted  current,  1, 
and  accelerating  voltage  as 


F  =  I 


(1-4) 


The  thrust  produced  is  directly  proportional  to  the  emitted  particle  current 
and  inversely  proportional  to  the  specific  charge.  Hence,  to  attain  8 
large  thrust  either  the  emitted  current  must  be  large  or  the  particle 
specific  charge  small  or  bo^h. 


B.  Power  Conversion  Efficiency 

An  electrostatic  thrustor  functions  by  converting  electrical  energy 

into  thrust.  The  efficiency  of  this  conversion  process  can  be  divided 

into  two  components.  The  first  concerns  the  transformation  of  electrical 

energy  into  beam  kinetic  energy,  and  the  second  the  relationship  between 

this  kinetic  energy  and  the  thrust  produced. 

The  formation  efficiency  of  an  electrostatic  thrustor  is  expressed 

in  terms  of  the  kinetic  energy  given  to  a  particle,  qVg ,  and  the  energy, 

W  ,  which  is  lost  in  forming  the  particle. 

L 


qV 

7)  = - 2— 

't  qV  +  W. 


i 
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The  formation  energy  consists  of  the  energy  needed  to  evaporate  the 
propellant,  control  the  propellant  flow,  support  the  ionization  itself, 
supply  heat  lost  through  conduction  and  radiation,  ohmic  losses,  and  space 
charge  effects.  The  formation  efficiency  is  plotted  3S  a  function  of  the 
accelerating  voltage  in  Fig.  (1-2).  In  order  to  maintain  a  high  formation 
efficiency  either  the  energy  loss  per  particle  must  be  low  or  the 
accelerating  voltage  high  or  both. 

The  thrust  producing  capability  of  a  charged  particle  beam  with  a 
specified  kinetic  energy  is  given  by  its  beam  conversion  efficiency.  The 
kinetic  power  of  the  beam  is  given  by 


P  =  1/2  mv2 


By  dividing  Eqn.  (1-3)  by  (1-6)  the  thrust  per  unit  power  is  obtained 

L-  -  ■"*  -  £  (1-7) 

P  1/2  mv*  ~  v  1  ; 

Most  physical  electrostatic  thrustors  do  not  expel  all  propellant 
particles  at  the  same  velocity.  Instead  the  exhaust  beams  possess  a 
velocity  distribution.  The  effect  of  a  velocity  distribution  is  to  reduce 
the  thrust  produced  per  unit  of  beam  power. 

If  the  system  has  a  velocity  distribution  given  by 


f(v)  =lFkvk 


where  F  is  the  fraction  of  particles  possessing  velocity  v  ,  then  the 

K  K 

thrust  per  unit  power  can  be  written  as 


IV 


(1-9) 


k  k  IVu 


f  gun  wffi l|  M  ifi'if  -*  :■  '■  "‘W 


FORMATION  EFFICIENCY 


6 


ACCELERATING  POTENTIAL  IN  KILOVOLTS 


FIG.  1-2  FORMATION  EFFICIENCY  AS  A  FUNCTION  OF  ACCELERATING 
POTENTIAL  WITH  FORMATION  ENERGY  AS  A  PARAMETER 
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For  a  distribution 
reduces  to  Eqn.  (1-7),  a 
multi-velocity  distribut 
given  by 


containing  only  c  single  velocity,  Eqn.  (1-9) 
nd  the  degradation  in  efficiency  produced  by  a 
ion,  or  beam  conversion  efficiency,  is  seen  to  be 


(ivj2 

XVk2 


(1-10) 


It  immediately  follows  that  the  summations  in  the  above  equation  can 
be  replaced  by  statistical  average  values 


(1-11) 


and  finally,  using  Eqn.  (1-2),  the  averaged  velocities  can  bo  replaced 
with  the  corresponding  averaged  specific  charges. 


(1-17) 


Accordingly,  in  an  electrostatic  thrustor,  the  beam  conversion 
efficiency  is  a  function  of  the  variance  of  the  specific  charge  distribu¬ 
tion  of  the  expelled  beam.  A  system  with  only  one  specific  charge  present 
produces  the  maximum  thrust  per  unit  of  beam  power;  as  the  distribution 
broadens  the  efficiency  is  reduced.  For  a  Maxwellian  distribution,  the 
beam  conversion  efficiency  is  approximately  85%. 

C .  Optimum  Exhaust  Velocity 

The  terminal  velocity  u  of  a  space  vehicle  in  a  drag  free,  gravity 
free  space  is  related  to  the  exhaust  velocity  of  .ts  propellant  v,  and 


8 


its  total  mass  ratio  M^/  by  Tsiolkovskii ' s  equation. 


(1) 


u  =  v  In  (  M  /M  ) 
o  e 


(1-13) 


where  M  ,  the  total  iritial  mass  and  M  ,  the  burnout  or  final  mass  are 
o  c 

given  by 


M  =  M  +  M  +  M, 
o  p  w  1 


M  =  M  +  M, 
e  w  1 


(1-14) 


where  M  is  the  propellant  mass,  M  the  power  plant  mass,  and  M,  the  pay- 

p  w  1 

load  mass . 

By  eliminating  the  power  plant  mass  from  Eqn .  (1-13),  the  pnyioad 

ratio  M  /  M  can  be  determined  as  a  function  of  the  exhaust  velocity.  The 
1  o 

propellant  flow  rate  m,  propulsion  time  tQ,  and  propellant  mass  ore 
related  by 


M  =  mt 
P  c 


(1-15) 


Hence,  the  beam  kinetic  power  can  be  written  as 

.2  M  v2 

P  _  in^_  =  _e — 

2  2t 


(1-16) 


The  beam  power  per  unit  of  power  plant  mass  is  defined  as  the  specific 
power  a  of  the  system 


a  = 


(1-17) 


By  combining  Eqns.  (1-13)  to  (1-17)  the  payload  ratio  M^/  is 
obtained  as  a  function  of  the  exhaust  velocity. 


M  2 

-rf-  =  EXP( -u/v)  -  -i— 
M  2Qrt 

o  o 


(1  -  EXP( -u/v) ) 


0-18) 


The  payload  ratio,  unlike  the  total  mass  ratio,  depends  not  only  on 
the  exhaust  and  terminal  velocities,  but  also  on  the  specific  power  ana 
total  propulsion  time. 

For  a  given  terminal  velocity,  specific  power,  and  propulsion  time 
the  payload  ratio  can  be  maximized  as  a  function  of  the  exhaust  velocity. 
The  location  of  the  muxlmum  is  found  bv  differentiating  Eqn .  (1-18)  with 
respect  to  the  exhaust  velocity.  The  optimum  exhaust  velocity  is  given  by 


2at 


2^  =  ^  <EXP(u/v)  -  1)  -  1 


(1-19) 


and  the  corresponding  maximum  payload  ratio  is 


M  1  -  (u/2v)  (1  +  v  /20Ct  ) 
_ i_  _ o_ 

Mo  1  +  (u/2v) (1  +  2 at  /v2) 

O 


(1-20) 


Consequently,  for  a  given  terminal  velocity  the  payload  ratio  can 
be  maximized,  likewise  at  a  specified  payload  ratio  the  terminal  velocity 
can  be  optimized.  The  payload  ratio  is  plotted  as  a  function  of  exhaust 
velocity  with  terminal  velocity  as  a  parameter  in  Fig.  (1-3).  All 
velocities  were  normalized  by  setting  2Ctt^  =  1. 


1 . 3  Comparison  of  Electrostatic  Propulsion  Systems 

In  the  previous  section  the  p.iysics  of  electrostatic  energy  conversion 
was  outlined.  Various  parameters  of  a  charged  particle  beam  were  deemed 
important  to  its  overall  performance  as  a  means  of  converting  electrical 
energy  into  thrust.  First,  the  formation  energy  of  a  particle  must  be 
small  compared  to  the  kinetic  energy  it  acquires,  second,  the  variance  of 
the  particle  specific  charge  distribution  should  be  small;  third,  the 


total  emitted  current  should  be  large;  and,  fourth,  the  exhaust  velocity 
should  be  optimized  for  cither  maximum  payload  ratio,  or  maximum  terminal 
velocity  or  other  mission  parameter. 

Electrostatic  thrusters  can  be  divided  into  two  classes  -  ionic  and 
heavy  particle  propellant.  Each  scheme  will  be  discussed  and  the  advan¬ 
tages  and  disadvantages  of  both  systems  will  be  compared  using  the  four 
criteria  listed  above. 


A .  Ionic  Propellant 

Ionic  propellant  electrostatic  thrustors  function  in  two  discrete 
steps.  First  an  ion  is  produced  by  one  of  several  known  ionization 
techniques,  and  then  it  is  accel.  rated  to  a  given  velocity  by  an  electro¬ 
static  field.  Many  different  methods  of  ionizing  atoms  have  been  developed 
since  Kutherf ord^  and  Thomson/5^  made  the  original  studies  on  the  con¬ 
duction  of  electricity  through  gases.  Langmuir"^  and  Taylor^7  performed 
much  of  the  original  experimental  and  theoretical  work  on  surface  ion 

(  g) 

sources.  In  19-10,  Finkelstein  examined  the  ionization  of  gases  by 
high  energy  electrons  and  developed  the  powerful  electron  bombardment 
ion  source.  A  third  type  of  ion  source  known  as  the  duoplasmatron  utilizes 
a  magnetic  pinch  effect  to  produce  a  very  dense  plasma;  it  was  conceived 

(9) 

by  von  Ardenne  in  1957. 

The  surface  ionization  source  was  the  first  of  the  above  to  be  pro- 
(10) 

posed  for  ion  engines  and  is  still  the  most  commonly  used  technique. 
Surface  ion  sources  utilize  the  fact  that  a  high  work  function  metallic 
surface  may  present  a  potential  well  to  the  electrons  of  a  low  ionization 
potential  atom.  Ionization  occurs  when  an  atom  strikes  the  metallic 


lurlacc.  This  surface  ionization  proceeds  in  two  discrete  steps.  First 


on  atom  approaching  the  surface  is  polarized  by  image  forces  between  its 
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nucleus  and  free  electrons  in  the  metal.  It  will  adhere  to  the  surface 
by  this  coulomb  attraction  os  an  atom.  When  the  temperature  of  the  sub¬ 
strate  is  hot  enough  the  atom  will  evaporate.  If  the  ionization  energy 
of  the  evaporating  atom  is  greater  than  the  work  function  of  the  substrate, 
the  atom  will  evaporate  as  a  neutral;  but,  if  the  work  function  is  greater 
than  the  ionization  energy  the  electron  is  trapped  by  the  substrate  and 
the  adsorbed  atom  evaporates  as  an  ion.  If  the  substrate  temperature  is 
below  a  critical  value,  the  atoms  will  evaporate  as  neutrals  and  above 
that  temperature  they  will  evaporate  as  ions. 

Much  work  has  been  spent  recently  on  the  investigation  of  the  ioni¬ 
zation  properties  of  cesium  on  porous  tungsten^1^  Cesium  is  the  heaviest 
element  that  is  useful  in  a  surface  ion  source  and  thus  is  the  most 
efficient  for  electric  thrustors.  Cesium  on  tungsten  ion  sources  have  a 
theoretical  ionization  efficiency  of  99.5%  and  laboratory  experiments 
have  yielded  efficiencies  in  this  range.  Hence,  such  a  source  produces 
over  99%  ions  and  less  than  1%  neutral  atoms,  and  the  variance  of  the 
corresponding  specific  charge  distribution  is  very  low. 

One  of  the  main  limitations  of  cesium  systems  is  the  energy  loss 

o 

per  ion  produced.  Surface  temperatures  in  the  1100  -  1500  K  range  are 
necessary  for  efficient  ion  production,  and  hence  heat  losses  by  conduc¬ 
tion  and  radiation  are  large  and  may  amount  to  more  than  1000  ev  per  ion. 
Therefore,  in  order  to  obtain  a  high  formation  efficiency  the  accelerating 
voltage  must  be  high,  but  this  in  turn  produces  a  high  exhaust  velocity. 
Thus,  for  many  missions,  a  cesium  ion  source  cannot  simultaneously  attain 
a  high  formation  efficiency  and  an  optimum  exhaust  velocity. 


B .  Heavy  Particle  Propellant 


The  basic  limitation  in  formation  efficiency  of  ionic  propellant 
systems  can  be  overcome  by  using  a  heavy  particle  propellant.  If  the 
maps  of  the  propellant  particles  is  increased  a  higher  accelerating  poten¬ 
tial  can  be  used  to  attain  a  desired  exhaust  velocity.  Hence,  by  varying 
the  specific  charge  of  the  propellant,  the  two  desirable  conditions  of 
a  high  accelerating  potential  and  an  optimum  exhaust  velocity  can  be 
achieved  simultaneously. 

The  most  widely  used  and  most  efficient  method  of  producing  heavy 
particles  is  by  the  eiectrohvdi odynamic  (EHD)  spraying  of  liquids  from 
electrified  capillaries.  By  placing  a  small  diameter  capillary  at  a  high 
potential  (either  positive  or  negative)  relative  to  a  nearby  grounded 
electrode  a  high  intensity  electric  field  can  be  produced  at  the  tip  of 
the  capillary.  If  a  liquid  is  forced  through  the  capillary  it  will  form 
a  hemispherical  meniscus  at  the  capillary  tip  and  a  surface  charge  will 
be  induced  on  it  by  the  electric  field.  If  the  outward  stress  of  the 
surface  charge  is  greater  than  the  surface  tension  forces  holding  the 
meniscus  together,  it  will  disrupt  and  a  spray  of  charged  particles  will 
be  produced. 

The  EHD  spraying  of  liquids  has  been  investigated  since  1745,  when 

(12) 

Bose  observed  the  formation  of  water  threads  at  the  tip  of  an  electri¬ 

fied  capillary.  The  physics  of  the  spraying  process  was  first  formulated 

(13) 

by  Rayleigh  in  1882  when  he  determined  the  stability  criterion  for  a 

charged  conducting  droplet.  Early  experimenters  believed  that  only  low 

conductivity  liquids  could  be  electrically  sprayed;  hence,  most  experiments 

(14) 

were  performed  with  insulating  liquids.  The  charge  and  size  of  indivi¬ 

dual  oil  droplets  electrically  sprayed  in  a  vacuum  were  measured  by 
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On  the  basis  ol  these  theoretical  studies  and  past  research  experience, 
studies  of  working  fluids  were  undertaken  in  an  attempt  to  optimize  the 
production  of  heavy  particles  and  eliminate  the  formation  of  ions.  In 
general,  the  EHD  spraying  of  non-conducting  liquids  has  produced  heavy 
particles  in  the  size  range  useful  for  the  electric  propulsion  systems. 
Unfortunately,  these  charged  particle  beams  have  only  been  produced  at  low 
current  levels  and  with  non-ideal  specific  charge  distributions.  Thus,  by 
using  heavy  particles  a  high  formation  efficiency  can  be  obtained,  but 
only  at  the  expense  of  a  low  emitted  current  and  a  less  than  100%  beam 
conversion  efficiency. 

1  .4  Proposed  Research 

Heavy  particle  systems  are  capable  of  high  formation  efficiencies 
and  ionic  systems  can  produce  high  beam  currents  at  high  beam  conversion 
efficiencies.  Heavy  particle  systems  were  first  introduced  as  a  means  of 
attaining  a  high  formation  efficiency  by  increasing  the  kinetic  energy 
given  to  the  propellant  particles  while  maintaining  a  near  optimum  ex¬ 
haust  velocity.  The  mere  direct  approach  of  decreasing  the  total  energy 
expended  in  producing  an  ion  has  met  with  little  success.  Almost  all  of 
this  energy  loss  in  surface  ionization  systems  can  be  attributed  to  heating 
effects  in  the  entire  system.  If  an  ionization  source  could  be  devised 
that  functions  at  low  temperatures,  the  formation  efficiency  could  be 
greatly  increased. 

One  method  of  producing  ions  at  room  temperature  is  by  the  field 
evaporation  of  a  surface.  A  sufficiently  strong  electric  field  applied 
to  a  metal  surface  in  such  a  direction  as  to  make  the  metal  positive 
produces  an  emission  of  positive  ions.  Field  evaporation  was  first 
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observed  by  Muller  In  1956^°^  while  experimenting  with  a  field  ion 
microscope  (FIM) .  Field  evaporation  techniques  have  been  utilized  as 
a  means  of  producing  very  smooth  emitters  for  the  FIM.  Since  any  irregu¬ 
larities  on  the  surface  lead  to  locally  high  electric  field  intensities, 
they  are  preferentially  evaporated  and  the  end  iorm  is  atomically  smooth. 

By  further  increasing  the  applied  electric  field  intensity  the  emitter 

can  be  disected  a  monolayer  at  a  time  and  hence  the  "inside"  of  the 

(21 ) 

emitter  can  be  examined. 

If  the  field  evaporated  ions  originate  in  an  adsorbed  layer  or  are 
produced  by  the  emitter  material  itself,  the  ion  source  will  be  weak  and 
rapidly  exhausted.  Muller  has  produced  steady  ion  currents  in  the 

(22 

microampere  range  by  supplying  the  ionizable  matter  in  the  gaseous  phase. 
The  ionization  of  gaseous  elements  at  an  emitter  held  at  a  high  potential 
is  the  basis  of  the  FIM. 

In  1967,  Swatik^  studied  the  EHD  spraying  of  the  eutectic  alloy 
of  gallium  and  indium  and  observed  that  under  the  appropriate  operating 
conditions  a  beam  consisting  totally  of  ions  could  be  produced  by  field 
evaporation  of  the  working  fluid.  Thus,  it  was  shown  for  the  first  time 
that  an  EHD  spraying  system  could  be  optimized,  not  for  the  production  of 
heavy  particles,  but  for  the  production  of  ions'.  In  such  a  system,  the 
ionizable  matter  supplied  in  the  liquid  phase  and  the  source  can  be 
operated  indefinitely  at  high  current  levels  because  the  evaporated  mass 
can  be  constantly  replenished. 

Thus,  it  has  been  the  goal  of  this  research  to  further  study  the  field 
evaporation  of  liquids  by  EHD  spraying  techniques.  Ideally  such  a  source 
of  charged  particles  would  have  a  very  low  energy  loss  per  particle 
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produced,  a  very  narrow  specific  charge  distribution,  and  a  high  beam 
current.  Hence,  it  would  possess  the  advantages  of  both  heavy  particle 
and  ionic  propellant  systems  and  would  represent  a  highly  efficient  means 
of  converting  electrical  energy  into  thrust. 

i 

The  wor  ing  fluid  used  for  field  evaporation  experiments  must  possess 
a  low  vapor  pressure  to  retard  evaporation  in  low  pressure  environments;  ■ 

a  high  surface  tension  to  maintain  stability  against  the  action  of  the 

f 

applied  electric  field;  and  a  high  conductivity  to  allow  for  the  production  j 

| 

i 

of  high  currents.  The  eutectic  alloy  of  76%  gallium  and  24%  indium  was 

chosen  as  the  working  fluid  for  this  study  on  the  basis  of  its  physical  I 

o  (24)  j 

parameters.  The  melting  point  of  the  alloy  is  15.7  C.  It  is  thus  a  j 

i 

liquid  at  room  temperature.  The  alloy  possesses  a  vapor  pressure  of  less 
-8  (25) 

than  10  torr  at  room  temperature.  The  surface  tension  and  conducti- 

(  26) 

vity  of  the  alloy  were  measured  by  Zrnic  and  SwatiK  and  are  628  dynes/ 

cm  and  3.4  x  104  mho/cm  respectively. 

The  first  phase  of  this  research  was  concerned  with  the  production 
of  a  stable  ion  emission  from  the  gal 1 ium- indium  alloy.  The  design  of 
the  capillary  and  accelerating  electrode  will  be  analyzed;  the  regulation 
of  the  accelerating  potential  will  be  determined  for  optimum  performance, 
and  the  stability  of  the  liquid  meniscus  will  be  examined. 

The  second  phase  was  the  analysis  of  the  particle  beam  with  respect 
to  its  performance  as  an  electrostatic  thrustor.  The  parameters  to  be 
measured  included  the  emitted  current  as  a  function  of  the  applied  electric 
field,  the  specific  charge  distribution  as  a  function  of  the  current,  and 
the  energy  loss  per  particle  as  a  function  of  the  current. 
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2 .  FIELD  EVAPORATION  OF  LIQUIDS 


The  mechanisms  involved  in  the  production  of  ions  from  an  electrified 
liquid  surface  will  be  discussed  in  this  chapter;  (1)  The  shape  which 
the  liquid  meniscus  should  take  under  the  influence  of  surface  tension 
and  the  applied  electric  field  will  be  determined;  (2)  An  approximate 
solution  of  Laplace's  equation  near  the  emitting  surface  will  be  formu¬ 
lated;  (3)  A  theory  explaining  field  evaporation  will  be  derived  from 
basic  principles;  and,  (4)  Several  possible  loss  mechanisms  in  the  beam 
formation  process  will  be  proposed. 


2.1  Shape  of  the  Meniscus 

The  action  of  electric  fields  on  liquid  surfaces  has  been  theoretically 
(13) 

studied  since  Rayleigh  derived  a  mathematical  stability  criterion  for 

electrified,  conducting,  liquid  spheres  in  terms  of  perturbations  on  their 

surfaces.  This  theory  explained  the  incipient  instability  of  a  drop  for 

infinitesimal  deformation  but  did  not  predict  or  explain  large  deformations. 

In  1964,  Taylor  studied  the  "Disintegration  of  Water  Drops  in  an  Electric 
••(27) 

Field"  and  found  that  a  conical  surface  could  exist  in  equilibrium 
under  the  influence  of  surface  tension  and  an  electric  stress,  and  thus 
approached  the  problem  of  large  scale  deformations. 

The  inward  stress  due  to  surface  tension  at  a  point  on  a  liquid 
surface  is  given  by 


r  -  y(L  + 1 ) 

\pi  p2/ 


(2-1) 


where^  is  the  surface  tension  constant  of  the  liquid  and  p^  and  p^  are 
the  two  principle  radii  of  curvature  at  the  point  in  question. 
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For  a  conical  surface  one  of  the  two  radii  is  infinite  at  all  points, 
and  Eqn .  (2-1)  becomes 


r  = 


r 


(2-2) 


It  is  evident  that  the  curvature  of  a  conical  surface  in  spherical 
coordinates  is  inversely  proportional  to  the  radial  coordinate  r.  The 
restoring  stress  of  surface  tension  is  then  also  inversely  proportional 
to  r,  and  an  outward  stress  normal  to  the  surface  that  can  be  in  equili¬ 
brium  with  the  surface  tension  must  also  be  inversely  proportional  to  r. 

The  stress  of  an  electric  field  normal  to  a  surface  is  a  function  of 

the  second  power  of  the  field  intensity.  This  Implies  that  the  potential 

1/2  -1 
must  vary  as  r  so  that  the  square  of  its  gradient  will  vary  as  r 

In  spherical  coordinates,  a  solution  to  Laplace's  equation  for  the 
region  bounded  by  a  conducting  conical  surface  and  one  of  its  equipoten- 
tials  is  given  by 


0 


+  B 


rrPn(cos6)j 


(2-3) 


where  P  (co  6)  is  the  Legendre  function  of  order  n.  This  geometry  is 
n 

shown  in  Fig.  (2-la). 

The  electric  field  that  is  required  for  the  equilibrium  of  the  conical 
surface  is  then  defined  by  the  potential 

1  /2 

0  =  Vo  +  B  r  P1/2(cos6)  (2_4) 


where  V  is  the  voltage  applied  between  the  conical  surface  and  the 
o 

equipotential . 

If  ©  =  0  is  the  conical  equipotential  surface  where  V  =  V  ,  then 
o  o 


1 


the  Lengendre  function  must  vanish  at  all  points  on  the  surface 


P1  /2^cos9^  =  ° 


(2-5) 


The  function  P  .  (cosB)  has  only  one  zero  in  the  range  0  <  0  <  «  at 

1  /  fc 

0  =  130.7°;  hence,  a  conical  rirface  can  exist  in  equilibrium  under  the 


influence  of  an  electric  field  if  it  has  a  vertical  angle  of  98.6 


2.2  Electric  Field  Calculations  for  Conical  Meniscus 

After  having  determined  that  a  conical  geometry  can  be  a  stable 
configuration,  it  is  necessary  to  analytically  determine  the  electric  field 
normal  to  the  surface.  The  electric  field  is  found  by  taking  the  negative 


gradient  of  the  potential  solution  given  in  Eqn.  (2-4)  and  applying  the 


appropriate  boundary  conditions 


P,  ,  (cos6)  0.5  a 


C  =  o  P 

V  1/2  1/2 

81 


,iM.*1W*"e)>  5, 

1/2  ^  (cos0) 


where  a,  is  the  interelectrode  spacing  at  9  -  0,  and  ^  and  £  are  unit 
vectors  in  the  r  and  9  directions  respectively.  At  r  =  0  this  solution 
gives  an  infinite  electric  field  intensity.  This  is  not  physically 
realizable.  To  make  the  problem  realistic,  the  con. cal  surface  can  be 
replaced  by  one  of  the  equipotent ia 1  surfaces  very  near  to  it,  defined  by 


r  =  a  (P,  .  (cos@) ) 
o  1/2 


(2-7) 


where  a  <<  a  .  The  equipotent  ia  1  approximation  is  shown  in  Fig.  (2-lb), 
o  1 

This  geometry  has  a  potential  solution  allowing  an  infinite  summation 
of  Legendre  functions,  but  the  approximation  is  m8de  that  the  only  term 


used  is  that  allowed  in  the  conical  solution  namely  n  =  1/2.  This  is 
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equivalent  to  assuming  that  the  equipotent ial  chosen  can  be  defined  by 
8  =  6  at  all  points  except  those  near  the  apex.  The  potential  is  then 
given  by 


(2-8) 


The  corresponding  electric  field  intensity  is 


Kence,  the  electric  field  intensity  at  the  tip  of  the  modified  conical 
surface  is 


£0.5  V 
(a  ,  o)  =  _ o_ 

°  1/2  1/2 


(2-10) 


a,  -  a 
1  o 


The  next  question  which  arises  concerns  the  relationship  between 

1/2  1/2 

the  ladius  of  curvature,  p,  of  the  equipotent iel  (r  vi/2  cose  =  80  ^ 

and  the  radial  distance  a  .  In  spherical  coordinates  the  curvature  of  a 

o 


surface  is  given  by 


(2-11) 


The  radius  of  curvature  at  a  point  is  defined  as  the  reciprocal  of  the 
curvature  at  that  point.  Evaluating  the  necessary  derivatives  on  the 
equipotent ial  surface  the  radius  of  curvature  at  the  tip  of  the  modified 
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conical  surface  is 


p(ao,n)  =  4aQ  (2-12) 

Using  Eqns.  (2-10)  and  (2-12)  the  electric  field  intensity  can  be 
obtained  as  a  function  of  the  radius  of  curvature  at  the  tip  of  the  modi¬ 
fied  conical  surface 


£<P0,0> 


0.5  V 

o 


(2-13) 


The  electric  field  intensity  is  plotted  for  various  values  of  applied 
voltage  at  an  assumed  interelectrode  spacing  of  0.25  cm  in  Fig.  (2-2). 


2 . 3  Field  Evaporation  Theory 

A  sufficiently  strong  electric  field  applied  to  a  metallic  surface  in 

such  a  direction  as  to  make  the  metal  positive  produces  an  emission  of 

positive  ions.  The  potential  energy  of  an  atom  or  ion  near  a  metal  sur- 

(23) 

face  can  be  approximated  by  a  single  curve.  If  at  a  given  distance  from 

the  surface  the  potential  energy  can  be  assumed  to  be  single  valued,  then 
energy  changes  only  occur  normal  to  the  surface  and  thermally  activated 
processes  can  be  explained  by  a  one-dimensional  energy  barrier.  Obviously 
such  a  model  is  valid  only  at  distances  which  are  large  compared  to  the 
dimensions  of  the  surface  structure. 

The  potential  energy  of  an  atom  as  a  function  of  distance  x  from  a 
surface  i c  shown  in  Fig.  (2-3a),  At  small  distances  this  curve  is  domi¬ 
nated  by  the  Pauli  exclusion  principle  and  at  large  distances  Van  der  Waal 
attractive  forces  determine  the  potential.  The  minimum  in  the  potential 
energy  curve  corresponds  to  the  sublimation  energy  A  of  the  atom.  The 


WITH  ACCELERATING  POTENTIAL  AS  A  PARAMETER 
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FIG.  2-3  ATOMIC  AND  IONIC  POTENTIAL  ENERGY  AS  A  FUNCTION 
OF  DISTANCE  FROM  A  METALLIC  SURFACE 
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corresponding  curve  for  an  ion  is  shown  in  Fig.  (2-3b)  .  The  energy  needed 

tc  ionize  an  atom  in  the  gas  phase  is  where  n  is  the  number  of 

"  n 

elemental  charges  of  the  ion  and  I  is  the  nth  ionization  potential  of  the 

n 

atom.  If  the  electrons  produced  by  the  ionization  are  returned  to  the  metal 
surface,  an  energy  n0Q  is  recovered.  The  quantity  is  the  absolute  work 

function  of  the  metal  which  is  the  energy  difference  between  an  electron 
at  infinity  and  an  electron  at  the  top  of  the  Fermi  surface  of  the  metal. 

The  zero  level  for  the  potential  energy  of  an  ion  must  therefore  be  at  a 
distance  ^1  -  n0Q  above  the  zero  level  for  the  corresponding  atom.  The 

potential  energy  curve  for  the  ion  is  determined  by  two  effects.  Again  at 
small  distances  the  Pauli  exclusion  principle  dominates,  but  at  large 
distances  the  ion  is  attracted  by  its  negative  image  in  the  metal  (Schottky 
image  potential)  and  a  minimum  in  the  potential  curve  is  produced. 

The  application  of  an  electric  field  to  the  specimen  alters  the 
potential  energy  curve  of  the  ion.  The  potential  energy  at  any  point,  x, 
is  reduced  by  n£ex  where  £  is  the  applied  field  and  e  is  the  charge  of 
an  electron.  The  effect  of  the  electric  field  is  to  produce  the  maximum 
in  the  ion  potential  energy  curve  shown  in  Fig.  (2-3c) .  At  a  distance  x 
from  the  metal,  the  ion  experiences  a  Schottky  image  force 


4ne  (2x) 2 
o 


(2-14) 


where  6  is  the  capacitivity  of  free  space, 
o 

By  including  the  effect  of  the  applied  field  and  assuming  that  the 
ion  is  sufficiently  far  from  the  metal  to  ignore  the  repulsion  potential, 
the  potential  energy  of  the  ion  is  given  by 


Uj  (  x) 
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(2-15) 
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The  maximum  of  this  potential  energy  curve  can  be  found  by  differen¬ 
tiating  Eqn .  (2-15)  with  respect  to  x.  The  position  of  the  maximum  is 


iound  to  occur  at 


tne 

is^rr 


(2-16) 


end  the  energy  at  this  maximum  is 


U.  (x  )  = 
I  max 


M 


(2-17) 


which  is  the  image  or  Schottky  hump  and  is  independent  of  x  as  long  as  the 
repulsion  potential  can  be  ignored. 

If  the  ionic  potential  curve  is  superimposed  on  the  atomic  potential 

(29) 

curve  several  possibilities  exist.  First,  if  the  ionization  potential 

is  sufficiently  low,  the  atom  may  be  adsorbed  at  the  surface  as  an  ion. 

The  ionic  potential  curve  lies  below  the  atomic  potential  curve,  and  field 
evaporation  takes  place  over  the  image  potential  hump  described  by 
Eqn.  (2-17)  and  shown  in  Fig.  (2-4a).  Second,  if  the  ionization  potential 
is  higher  the  atom  is  adsorbed  as  a  neutral  and  is  ionized  when  the  ionic 
and  atomic  curves  intersect.  If  this  intersection  is  to  the  left  of  the 
image  hump  as  shown  in  Fig.  (2-4b) ,  the  image  potential  theory  is  valid 
and  for  either  of  the  above  two  cases  the  activation  energy  for  field 


evaporation  Q  is  given  by 


Q  =  Q  - 
o 


-|] 

n3e3£ 

~4r,e 

l  °  J 


(2-18) 


where  Q  =  A  +  I  -  n  0 
o  1  *■  "  n  o 

Third,  at  still  higher  ionization  potentials  the  atomic  and  ionic 
curves  may  intersect  to  the  right  of  the  image-potential  hump,  as  shown 


in  Fig.  (2-4e),  and  the  image-potential  theory  is  no  longer  vslid.  In 


KgMtibHMlft . fcuoaii—i 


POTENTIAL  ENERGY  CURVES 
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k  is  Boltzmann’s  constant,  and  T  Is  the  temperature.  Substituting 
Eqn .  (2-18)  into  the  above,  the  time  constant  for  evaporation  under  the 


Influence  of  an  electric  field  is  then  given  by 


(2-22) 


The  total  current  emitted  by  field  evaporation  is  equal  to  the  current 
per  source  multiplied  by  the  number  N  of  lattice  sites  undergoing  evapora¬ 
tion  . 

1  «  Ne  (2-23) 

t 

From  Eqns  -  (2-22)  and  C~23)  the  emitted  current  can  be  found  ss  a  function 
of  the  applied  electric  field. 


I 


(2-24) 


As  the  electric  field  intensity  is  increased,  the  current  increases  due 
to  two  simultaneous  effects.  First,  the  time  constant  for  evaporation 
decreases;  and  second,  the  number  of  sites  contributing  to  the  total 
current  output  increases. 
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2 .4  Energy  Loss  In  Ion  Beam  Formation 

Many  possible  loss  mechanisms  exist  in  the  production  of  charged 
particle  beams  by  field  evaporation.  Three  sources  of  loss  will  be 
examined.  (1)  The  binding  energy  of  the  evaporated  atoms;  (2)  Ohmic 
losses  in  the  liquid  meniscus;  and  (3)  Space  charge  effects. 

A .  Binding  Energy 

The  energy  necessary  to  field  evaporate  an  atom  as  given  hv 
Eqn.  (2-18)  is  composed  of  three  terms.  The  height  of  this  energy  hump 
is  equal  to  the  sublimation  energy  of  the  atom,  plus  its  ionization 
energy,  minus  the  electron  activation  energy.  Therefore,  the  ions  are 
not  accelerated  through  the  entire  potential  applied  to  the  surface. 

A  fraction  of  this  energy  is  required  to  produce  the  particles.  The 
potential  that  is  needed  to  produce  an  ion  is  given  by 

V,  =  eQ  (2-25) 

f  o 

where  is  the  energy  required  to  create  an  ion.  For  both  Gallium  and 

Indium  Q  is  less  than  5  ev;  hence,  the  potential  drop  V.  is  insignificant 
o  i 

compared  to  the  kinetic  energy  that  the  ions  acquire. 

B.  Ohmic  Losses 

The  finite  conductivity  of  the  working  fluid  produces  a  non-zero 
value  of  resistance  in  the  spraying  meniscus.  The  physical  conical 
meniscus  has  a  finite  radius  apex. 

This  geometry  can  be  approxim^ ced  by  a  truncated  cone.  The  resis¬ 
tance  of  such  a  conical  volume  is  given  by 


32 


R,  p 


r^/tana 


dz 


itzs(tana)4 


(2-26) 


'r^/tana 


where  p  is  the  resistivity  of  the  working  fluid,  a  is  the  semivertical 
angle  of  the  conical  surface,  z  is  a  "dummy"  variable  of  integration,  and 
and  are  the  minimum  and  maximum  radii  cf  the  conical  surface 
respectively. 

Equation  (2-26)  is  easily  integrated  and  assuming  ^<<^1  t^le 
following  expression  is  obtained  for  the  resistance  of  the  truncated 
cone . 


R  ^ 


r^  Mtana) 


(2-27) 


For  the  Gallium-Indium  eutectic  working  fluid  and  a  semivertical  angle 


of  49.3  ,  Eqn.  (2-27)  becomes 


R  = 


8.1x10 


-6 


(2-28) 


The  potential  drop  in  the  liquid  is  then  given  as  the  product  of 
the  emitted  current  and  the  resistance  of  the  meniscus 


V  =  1R 

7. 


(2-29) 


,-8 


For  a  current  of  500 ymA  and  mini  urn  radius  of  x0  meter,  the  resulting 
potential  drop  is  le:s  than  one  volt.  Consequently,  for  the  Gallium- 
Indium  working  fluid  ohmic  losses  are  negligible. 
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C .  Space  Charge  Effects 

The  presence  of  positive  charges  in  the  vicinity  of  the  field  evapor¬ 
ating  surface  alters  the  potential  distribution  of  the  applied  field. 

Since  r>  cnrrge  placed  in  a  volume  already  containing  a  net  charge  density 
possesses  a  well  defined  potential  energy,  it  is  necessary  for  the  emitted 
particles  to  convert  some  of  their  available  kinetic  energy  into  potential 
energy  in  order  to  be  placed  into  the  beam.  To  find  the  magnitude  of  this 
potential  barrier,  it  is  assumed  that  the  emitted  positive  particles  form 
a  cylindrical  beam  with  a  constant  current  density  J  over  its  cross- 
sectional  erea. 

The  current  density  J,  volume  charge  density  p ,  and  particle 
velocity  v  are  related  by 


J  = 


pv 


(2-30) 


The  particle  velocity  in  turn  is  related  to  the  accelerating  potential 
by  Eqn.  (1-2).  Poisson's  equation  in  cylindrical  coordinates  can  be 
written  in  the  form 


d2V  +  1  dV 


dr 


r  dr 


Jr  (2q  V>-‘/2 

c  s 

o 


(2-31) 


This  equation  can  be  solved  more  easily  by  normalizing  the  variables  V 
(33) 


and  r  . 


Let  the  potential  at  the  center  of  the  beam  be  V  ,  and  the 


potential  at  any  other  radius  be  V  +  V(r) .  The  potential  may  be  written 
in  a  dimensionless  form 


V_ 

V 


1  +  w 


(2-32) 
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where  W  =  v(r)/v  . 

c 

Using  the  normalized  potential,  Poisson's  equation  becomes 


dfw  1^  dW 
dr2  +  r  dr 


J(2q  ) 


-1/2 


e  v  # 

o  c 


a  +  w) 


-1/2 


(2-33) 


The  radial  coordinate  may  also  be  put  in  a  dimensionless  form  by 
making  the  cori espondence  R  —  r/g,  where  g  is  suitable  normalizing  radius 
defined  by 


2 

g 


e  V  1/2(2q  )1/2 

_ o  c _ JS _ 


(2-34) 


Hence,  the  transformed  Poisson  equation  can  be  written  as 


d2W  j  1  dW 
di5  +  R  dR 


(1  +  W) 


-1/2 


(2-35) 


Since  there  is  no  simple  closed  i orm  solution  to  Eqn.  (2-35),  it  is 
solved  most  easily  by  a  power  series  expansion.  Noting  that  W  is  zero  at 
R  =  0,  and  that  W  is  symetnc  with  respect  to  R,  one  can  assume  a  solu¬ 
tion  of  the  form 


W 


V 


V 


V 


(2-36) 


Substituting  the  above  expansion  into  Eqn.  (2-35)  and  expanding 
the  right  hand  side,  the  solution  obtained  is 


2  2 

R  .  R 
-  (1  4 - + - +  .  .  . 

4  32  2304 


(2-37) 


For  R’<C  0.1  the  first  term  of  this  series  approximates  the  true 
solution  to  an  accuracy  of  greater  than  99%.  The  increase  m  potential, 
due  to  space  charge  at  the  center  of  the  positive  ion  beam,  relative  to 
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the  perimeter  of  the  beam  is  given  by 


V 

B 


_ I _ 

4e  (2q  V  )1/2 

o  sc 


(2-38) 


Hence,  to  first  order,  the  effect  of  space  charge  is  independent  of 
the  beam  radius  and  directly  prop  irtional  to  the  total  current.  For  a 
current  of  500  mA,  a  voltage  of  10  kv,  and  a  specific  charge  of  5  x  10s 
coul/Kg,  the  potential  due  to  the  space  charge  is  approximately  150  volts. 
This  value  is  1.5%  of  the  total  assumed  potential  of  10  kv. 

Therefore,  in  considering  the  three  loss  mechanisms  discussed  in 
this  section  it  is  obvious  that  only  space  charge  effects  would  produce 
measureable  energy  losses.  The  effects  of  the  space  charge  would  be  to 
produce  a  transverse  acceleration  of  the  ions  in  the  beam. 


EXPERIMENTAL  TECHNIQUES 


The  experimental  apparatus  and  instrumentation  will  be  described 
and  relevant  design  criteria  expounded  upon.  The  discussion  will  be 
subdivided  into  three  categories:  (1)  The  vacuum  system;  (2)  The  EHD 
sDraying  apparatus;  and  (3)  The  particle  detection  instrumentation. 

3 . 1  Vacuum  System 

All  experimentation  was  performed  in  an  ambient  pressure  of  less 

_  g 

than  10  torr  that  was  maintained  by  an  oil  diffusion  pump  and  a  mecha¬ 
nical  roughing  and  fore  pump. 

The  vacuum  chamber  consisted  of  a  custom  made  six  inch  industriaj 
glass  pipe*  tee  with  an  additional  four  inch  port  mounted  perpendicular 
to  the  plane  of  the  tee.  The  three  six  inch  ports  available  were  used 
for  mounting  the  spraying  apparatus,  mounting  the  detection  apparatus, 
and  pumping  respectively.  The  four  inch  port  was  utilized  for  an  optical 
monitoring  system.  A  photograph  with  an  accompanying  schematic  diagram 
of  the  experimental  set-up  are  shown  in  Fig.  (3-1). 

3 .2  Spraying  Apparatus 

The  beam  formation  equipment  consisted  of:  (1)  A  set  of  electrodes 
composed  of  the  spraying  capillary  and  accelerating  electrode,  (2)  A 
positive  displacement  feed  system  to  control  the  flow  of  working  fluid 
through  the  capillary;  and  (3)  The  accelerating  voltage  supplies.  The 
spraying  apparatus  is  shown  in  Fig.  (3-2). 

*  Pyrex  "Double  Tuff"  Brand,  Corning  Glass  Works,  Corning,  New  York 
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A .  High  Voltage  Electrodes 

The  electrode  geometry  consisted  of  the  capillary  (anode)  and  a 
concentric,  Pi  anner  accelerating  electrode  (cathode)  through  which  the 
capillary  protruded.  The  design  of  the  capillary  was  paramount  to  reali¬ 
zation  of  the  proper  mode  of  spraying  (i.e.  field  evaporation  mode). 
Relevant  design  factors  included  the  material,  tip  finish,  and  size  of 
the  capillary. 

The  choice  of  capillary  material  was  based  on  its  susceptibility 

to  wetting  by  the  ga 1 1 ium- l nd ium  eutectic.  Traditionally,  stainless 

steel  capillaries  have  been  employed  in  most  EHD  spraying  experiments 

because  of  their  ease  of  fabrication  and  accessibility  (standard  medical 

type  hypodermic  needle  stock).  1’nfortunatelv ,  the  gallium-indium  alloy 

does  not  properly  wet  stainless  steel.  Without  sufficient  wetting,  the 

predicted  conical  meniscus  is  not  produced.  Instead,  as  the  applied 

electric  field  is  increased,  the  entire  meniscus  is  pulled  off  of  the 

capillary  in  a  large  "blob".  The  desired  wetting  was  achieved  by  using 
(32) 

nickel  capillaries  that  were  fabricated  from  seamless  (250  OD 

x  125  ID)  nickel  tubing.* 

The  tip  of  the  capillary  har'  have  a  smooth  highly  polished 
surface.  A  small  burr  on  the  tip  could  lead  to  a  locally  high  electric 
field  and  perturb  the  spray.  The  raw  nickel  tubing  was  cut  to  approxi¬ 
mately  the  desired  length  with  diagonal  cutters.  This  procedure  crimped 
the  ends  of  the  tubing  shut.  They  were  then  placed  in  a  high  speed 
jewelers  lathe.  By  applying  a  slight  lateral  pressure  to  the  rotating 
tubing,  a  fatigue  crack  was  produced  that  eventually  ruptured  the  tubing 


*  .Nickel  200  Brand  Tubing,  Superior  Tube  Co.,  Norristown,  Pennsylvania 
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leaving  an  open  end.  The  spraying  tip  was  then  polished  with  jewelers 
rouge.  A  photomicrograph  of  a  finished  capillary  is  shown  in  Fig.  (3-3). 

The  size  of  the  capillary  influenced  the  stability  of  the  meniscus. 
As  the  capillary  size  was  decreased  the  meniscus  became  more  stable.  The 
lower  limit  on  the  capillary  size  was  dictated  by  the  availability  of 
small  diameter  nickel  tubing,  the  ease  of  fabrication  of  the  finished 
capillary,  and  its  resistance  to  plugging. 

B.  Feed  System 

The  finished  capillaries  were  cemented  into  the  cannula  of  a  hypo¬ 
dermic  needle.  This  facilitated  attachment  to  a  syringe  that  functioned 
as  a  reservoir  for  the  working  fluid.  The  flew  of  liquid  through  the 
capillary  was  controlled  by  a  positive  displacement  feed  system.  The 
principle  components  of  this  system  were  a  guide  mechanism  and  a  lead 
screw  that  transformed  an  angular  motion  into  a  linear  displacement.  The 
lead  screw  was  rotated  by  a  chart  recorder  drive  at  speeds  as  low  as  2 
revolutions  per  day.  The  linear  motion  produced  was  transferred  to  the 
plunger  of  the  syringe  and  the  resulting  mass  flow  rate  could  be  varied 
from  10  5  to  10  kg/sec . 

C .  Supply  Volt  ges 

The  electric  field  normal  to  the  conical  liquid  meniscus  must  be  of 
sufficient  amplitude  to  support  the  conical  geometry,  but  not  so  great 
that  it  pulls  off  the  entire  meniscus.  Because  the  gallium-indium 
eutectic  is  a  very  good  conductor,  the  meniscus  forms  the  anode  of  the 
spraying  field  geometry,  and  any  motion  of  the  liquid  changes  the 
"constant"  in  the  voltage-f if  1 d  relationship.  From  Eqn .  (2-24)  it  is 
evident  that  the  spraying  current  varies  greatly  with  a  small  change  in 
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fiold  strength.  Hence,  it  is  necessary  to  have  .'tdbeck  in  the  accel¬ 
erating  voltage  supply  to  maintain  a  constant  field  (i.e.  constant  current) 
independent  of  changes  in  the  electrode  geometry. 

The  accelerating  voltage  was  obtained  from  a  regulated  30  KVDC 
supply.  The  circuit  shown  in  Fig.  (3-4)  was  used  to  regulate  the  emission 
current  level  by  controlling  the  bias  on  a  series  regulator  tube.  A  three 
digit-digital  voltmeter  (DVM)  was  used  to  monitor  the  accelerating  voltage. 
The  input  impedance  of  the  metering  circuitry  was  approximately  1  gigaohm 
and  thus  did  not  effect  the  current  regulation. 

The  field  evaporated  ions  upon  impact  with  metallic  parts  in  the 
vacuum  system  can  produce  secondary  electrons.  These  secondaries  can 
then  be  accelerated  to  the  capillary  (at  a  high  positive  DC  potential) 
and  heat,  it  and  even  melt  it  shut.  By  placing  the  accelerating  electrode 
at  a  negative  bias  of  300  volts,  any  secondaries  produced  "see"  this 
negative  potential  instead  of  the  positive  capillary  potential.  (This 
is  because  the  bias  field  is  fairly  uniform  and  the  capillary  field 
decreases  rapidly  with  distance.)  Consequently,  the  secondaries  are 
suppressed,  the  capillary  is  not  heated,  and  the  measured  capillary 
current  is  only  the  ion  current . 


3 . 3  Detection  Apparatus 

The  following  parameters  of  the  field  evaporated  \  icles  were  to 
be  measured:  (1)  The  beam  current  as  a  function  of  the  applied  field; 
(2)  The  specific  charge  distribution  as  a  function  of  the  beam  current, 
and  (3)  The  formation  energy  of  the  particles  as  a  function  of  the  beam 


current 
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FIG.  3-4  HIGH  VOLTAGE  CURRENT  REGULATOR 
AND  PULSER  CIRCUITRY 
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The  measurement  of  the  specific  charge  distribution  was  accomplished 
by  utilizing  two  specific  charge  selective  instiuments.  The  first,  a 
time-of-f light  mass  spectrometer  (TOFMS)  was  used  because  it  can  detect 
a  large  cross-section  of  the  ion  beam.  Unfortunately,  it  has  low 
resolution.  The  second  device  used  was  a  quadrupole  mass  filter.  It 
was  placed  in  tandem  with  the  TOFMS  and  its  operation  complemented  that 
of  the  TOF  device.  The  quadrupole  detects  only  those  ions  very  near  its 
axis  of  symmetry,  but  it  is  capable  of  high  resolution.  Thus,  the 
simultaneous  use  of  both  instruments  provided  an  overall  picture  of  the 
specific  charge  distribution.  Since  both  instruments  sense  the  specific 
charge  of  the  particles  and  not  their  mass,  the  mass  spectrum  was  found 
by  assuming  that  all  particles  carried  only  one  elementary  charge. 


2 


(3-1) 


where  is  a  geometric  constant  of  tne  device.  In  general,  the  voltages 
applied  to  the  electrodes  can  be  written  as 


V(t)  =  t  (VJ  +  V  cosvJt ) 
dc  ac 


and  the  corresponding  potential  distribution  ; n  the  xy  plane  is 


V2  — 

c,v)  =  (V  +  V  co  scat )  : - - 

dc  ac  £ 
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The  electric  field  intensity  between  the  electrodes  is  then 

r~  -i 

xa  -  y£ 

F  =  2(V  +  V  cosurc )  — - - - - ^  (3-3 
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o 

The  equations  of  motion  of  positive  ions  under  the  influence  of  the  above 
electric  field  are 
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By  making  the  correspondences 
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the  equations  of  motion  are  transformed  into  standard  Mathieu's  equations 


+  (a  +  2qcos29)  x  =  0 


-  (a  +  2qcos2S)  y  =  0 

aw 


The  general  solution  of  these  equations  takes  an  exponential  form  and 


there  exists  finite  limits  on  the  ranges  of  a  and  q  for  stable  trajectories. 
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By  varying  the  ratio  a/q  the  resolution  of  the  nass  filter  is  changed.  At 

a  =  .236  and  q  =  .706  only  one  mass  number  is  stable  and  maximum  resolu- 
(34) 

tion  is  obtained.  At  maximum  resolution,  the  stable  mass  number  A,  BF 

amplitude  V  and  RF  frequency  ca  sre  related  by 

V  =  7.21  f2(MHz)  r  2(cu)  A  (3-8) 

0  c  o 

Therefore,  by  measuring  the  amplitude  of  the  RF  voltage  and  its  frequency 
at  a  particular  operating  point,  the  stable  mass  number  is  determined. 

The  quadrupole  mass  filter  was  fabricated  from  0.250  inch  diameter 
centerless  ground  stainless  steel  rods  40  cm  long.  They  were  end  supported 
by  0.125  inch  diameter  synthetic  sapphire  ball  bearings.  This  support 
system  allowed  for  precise  alignment  of  the  rods  and  provided  electrical 
isolation.  The  entire  assembly  was  placed  in  a  conducting  cylinder  to 
shield  the  RF  voltages. 

The  power  supply  for  the  quadrupole  is  shown  in  Fig.  (3-6).  A  basic 
Ecceles  Jordan  configuration  was  used  with  a  tuned  circuit  plate  load  to 
produce  a  double  ended  oscillator,  operating  at  approximately  2  MHz.  The 
RF  signal  obtained  was  stepped  up  to  3  KV  peak  amplitude  by  an  air  core 
inductor.  The  DC  bias  for  the  quadrupole  was  obtained  by  rectifying  a 
part  of  the  RF  signal  and  then  adding  it  back  to  the  RF  in  the  output 
stage . 

The  stable  mass  number  of  the  quadrupole  can  be  varied  by  varying 
either  the  RF  amplitude  or  RF  frequency.  Amplitude  sweeping  was  chosen 
because  of  its  simplicity.  A  block  diagram  of  the  quadrupole  instrumenta¬ 
tion  is  shown  in  Fig.  (3-7).  The  plate  voltage  for  the  quadrupole  supply 
was  obtained  from  a  programmable  power  supply.  By  applying  a  ramp  input 
to  the  supply,  the  plate  voltage  could  be  swept  and  hence  the  RF  output 


QUADRUPOLE 


FIG.  3-7  BLOCK  DIAGRAM  OF  QUADRUPOLE 
MASS  FILTER  INSTRUMENTATION 
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could  be  varied  to  cover  a  range  of  stable  mass  numbers  ranging  1 rom  l  to 
420.  The  ion  current  passing  through  the  quadrupole  was  amplified  by  a 
high  gain  electrometer  circuit  and  the  output  of  the  electrometer  was 
used  to  drive  a  strip  chert  recorder. 
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The  current  density  as  a  function  ol  the  specific  charge  (specific 
charge  distribution)  is  related  to  the  current  densitv  as  a  function  of 
drift  time  by 


(3-12) 


Accordingly,  the  specific  charge  distribution  can  be  obtained  by  differ¬ 
entiating  the  beam  current  waveform. 

A  schematic  diagram  of  the  TOFMS  is  shown  in  Fig.  (3-b) .  The 

collector  was  formed  with  a  spherical  surface  with  a  radius  of  25  cm 

corresponding  to  the  collector-capillary  separation.  This  arrangement 

o 

provides  analysis  oi  all  ions  sprayed  into  a  cone  of  approximately  30 
apex  angle.  Ions  sprayed  at  a  greater  angle  hit  the  surrcunding  shield. 
These  ions  could  suffer  elastic  collisions  with  the  shield  and  still  hit 
the  collector,  or  they  could  produce  secondary  electrons  upon  impact.  To 
shield  the  collector  from  rebounding  ions  and/or  secondary  electrons  a 
shield  is  placed  between  the  capillary  and  the  collector  with  an  aperature 
of  proper  size  so  that  all  primary  particles  outside  the  cone  of  detection 
are  blocked  from  the  collector. 

The  third  element  of  the  TOFMS  is  an  electrically  floating  fine  wire 
grid  placed  near  the  collector,  that  shields  it  from  the  capillary.  By 
biasing  this  grid  negative  it  also  functions  to  suppress  secondary 
electrons  produced  by  ion  bombardment  of  the  collector  surface. 

The  ion  beam  was  pulsed  off  by  shunting  the  accelerating  voltage  to 
ground  potential  with  a  high  voltage  thyratron.  The  pulser  network  is 
shown  in  Fig.  (3-4).  The  fall  time  of  the  resulting  pulse  was  typically 
less  than  0.1  microsecond  and  was  much  less  than  the  average  ion  drift 
time  of  2.5  microseconds.  The  length  of  the  pulse  was  also  important. 


BEAM  COLLECTOR 

EXTERNAL  SHIELD  DIVIDER  GRID 


r>3 

Since  the  iced  system  ran  continuously,  a  long  accelerating  voltage 
pulse-off  could  produce  an  excess  ol  mass  at  the  capillary  and  disrupt 
the  spray.  This  problem  was  circumvented  by  making  pulses  of  relatively 
short  duration  typically  less  thun  10  ms.  This  was  accomplished  de¬ 
limiting  the  current  through  the  thyratron  to  a  value  less  than  its  hold¬ 
ing  current.  In  this  mode  of  operation,  the  thyratron  will  remain  "on" 
only  as  long  as  its  grid  is  biased  positive.  Thus  by  pulsing  the  grid, 
the  "on"  time  could  be  controlled. 

The  collector  current  was  monitored  by  a  cathode  ray  oscilloscope 
with  a  33  MHz  bandwidth.  The  capacitance  of  the  collector,  cabling,  and 
oscilloscope  input  was  approximately  100  pf  .  To  utilize  the  fast  rise 
time  of  the  oscilloscope  the  collector  was  shunted  to  ground  with  a 
1000  ohm  resistor  to  keep  the  time  constant  of  the  circuit  as  low  as 
possible . 

Data  was  taken  by  pulsing  the  accelerating  voltage  and  simultaneously 
triggering  the  esc i 1 1  /scope .  The  decay  of  the  collector  current  as  a 
function  of  time  was  recorded  by  an  oscilloscope  camera.  The  recorded 
current  curve  was  numerically  differentiated  to  obtain  the  specific 
charge  distribution. 

B.  Formation  Energy  Measurement 

The  energy  distribution  of  the  field  evaporated  ions  was  measured 
with  a  retarding  potential  analyzer.  The  basic  principle  of  the  cell  is 
to  produce  a  potential  barrier  of  sufficient  height  so  that  ions  with 
energy  lower  than  the  barrier  are  deflected  and  ions  with  higher  energies 
pass  unaffected. 

The  physical  retarding  potential  cell  utilized  the  TOFMS  drift  space. 
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An  additional  grid  was  inserted  between  the  capillary  and  the  collector 
on  to  which  the  retarding  potential  could  be  applied,  as  shown  in 
Fig.  (3-8) .  Tlie  retarding  potential  was  supplied  by  a  highly  regulated 
DC  supply  end  measured  to  an  accuracy  of  1  volt.  The  capillary  voltage 
was  measured  with  a  digital  voltmeter  that  was  calibrated  by  the  retarding 
potential  supply  to  an  accuracy  of  10  volts. 

If  the  energy  distribution  of  the  field  evaporated  particles  is 
given  by  the  function  F(E),  the  current  to  the  collector  is  equal  to 

c 

I  =  \  F(E)  dE  (3-13) 

where  V  is  the  retarding  potential. 

R 

Therefore,  the  energy  distribution  is  given  by  the  derivative  of 


the  collector  current. 


to  obtain  the  energy  distribution  function. 


•1  .  EXPERIMENTAL  RESULTS 


Various  mcasurwi’iits  were  made  of  the  parameters  of  the  field 
evaporation  mode  of  EHD  spraying.  These  data  were  taken  lor  two  purposes; 
(1)  to  corroborate  the  theory  explaining  this  physical  phenomenon  proposed 
in  Chapter  2;  and  (2)  to  determine  the  overall  efficiency  and  usefulness 
of  this  technique  us  it  .source'  of  charged  particles  for  ml crot hr ust or 
ut 1 1  l  za  t ion . 

Photomicrographs  ol  the  capillary  and  liquid  meniscus  were  made 
under  all  operating  conditions  to  aid  in  analyzing  the  operation  of  the 
ion  source,  and  to  determine  the  stability  of  the  Taylor  cone.  The 
current-voltage  characteristic  of  the  emitted  Ions  was  determined  and 
used  to  support  the  field  evaporation  theory  set  forth.  The  formation 
energy  of  the  particles  was  measured  at  various  current  levels  in  order 
to  determine  the  loss  mechanisms  in  the  particle  production. 

The  usefulness  of  the  field  evaporation  source  as  j  microthruster 
was  determined  by  analyzing  the  specific  charge  distribution  of  the  emitted 
particles  as  a  function  of  the  capillary  current.  Using  this  data  and 
the  current-voltage  c-ha  roe  t  er  l  si  i  e ,  the  thrust  and  beam  conversion 
efficient’.  i  the  source  could  be  determined.  Likewise,  the  formation 
efficiency  could  be  deduced  from  the  retarding  potential  data. 

4.1  S tabilization  of  the  Taylor  Cone 

The  field  evaporation  mode  of  EHD  spraying  was  experimentally  observed 
to  occur  whenever  the  Taylor  cone  was  produced.  The  formation  of  the  cone 
was  always  accompanied  by  an  emission  ol  ions  lion'  its  apex.  Unstable 
operating  points  resulting  in  oscillations  ol  the  conical  meniscus  produced 
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large  lowly  charged  particles.  The  cone  could  be  easily  reproduced  and 
was  indefinitely  stable  as  long  as  the  design  criteria  outlined  in 
Chapter  3  were  maintained.  The  two  most  important  parameters  affecting 
the  formation  of  the  cone  were  the  wetting  of  the  capillary  and  the  feed¬ 
back  control  (current  regulation)  of  the  accelerating  voltage. 

The  nickel  capillaries  promoted  wetting  by  the  ga 1 1 i um- i ndium  alloy 
and  were  pre-wet  v  - h  the  alloy  prior  to  insertion  in  the  vacuum  s.,  stem. 

Once  proper  wetting  was  achieved  the  capillary  remained  wetted  unless  the 
field  was  allowed  to  evaporate  off  all  of  the  working  fluid. 

The  effect  of  wetting  can  be  seen  in  the  sequence  of  photomicrographs 
shown  in  Fig.  (4-1).  A  250  j*m  O.D.  nickel  capillary  is  shown  in  Fig.  (4-la). 
In  Fig.  (4-lb)  a  sphere  of  working  fluid  has  been  created  at  the  tip  of 
the  capillary  by  the  feed  system.  It  does  not.  wet  the  capillary  and  thus 
does  not  flow  back  around  it.  The  applied  electric  field  distorts  the 
sphere  and  tries  to  pinch  it  off.  This  pinching  effect  reduces  the  mini- 
.•.-■i.  circumference  of  the  liquid  meniscus  and  hence  reduces  the  surface 
tension  forces  at  that  point  allowing  for  further  pinching  and  the  eventual 
breakoff  of  a  large  lowly  charged  droplet. 

Microphotographs  of  the  same  capillary  after  induced  wetting  are 
shown  in  Fig.  (4-lc).  In  this  case  the  surf  Bee  tension  forces  stabilizing 
the  meniscus  act  between  the  capillary  and  the  working  fluid.  As  the 
field  distorts  the  meniscus  the  surface  tension  forces  remain  constant 
and  support  the  field.  Further  increases  in  the  field  strength  can  then 
produce  the  Taylor  cone. 

The  regulation  of  the  accelerating  voltage  had  a  pronounced  effect 
on  the  stability  of  the  meniscus  The  system  was  started  up  by  operating 


58 


the  feed  system  until  a  hemispherical  meniscus  appeared  at  the  capillary 
tip.  This  meniscus  usually  had  a  larger  volume  than  that  needed  for  the 
stable  conical  geometry  and,  therefore,  the  meniscus  would  become  unstable 
and  large  particles  with  low  specific  charge  would  be  ejected  until  the 
proper  volume  of  working  fluid  remained.  In  this  mode  of  operation,  the 
application  of  an  unregulated  accelerating  voltage  would  result  in  the 
disruption  of  the  entire  meniscus.  By  use  of  current  regulation  the 
instability  could  be  controlled  until  a  stable  emission  was  obtained.  In 
this  mode  of  operation,  the  accelerating  voltage  was  continuously  varying 
over  a  dynamic  range  of  approximately  2  KV . 

By  changing  the  feedback  resistance  in  the  current  regulator  circuit 
the  capillary  current  could  be  continuously  varied  from  20  to  500  microamps. 
Once  the  conical  meniscus  was  stabilized,  the  applied  voltage  would  remain 
constant  to  within  100  volts  at  a  constant  current  level.  Without  the 
regulation  it  was  necessary  to  compensate  manually  for  the  si igntest 
movement  in  the  meniscus  in  order  tc  keep  a  steady  current  level.  Thus, 
the  current  regulation  was  found  to  aid  in  stabilizing  the  meniscus  and 
to  increase  the  dynamically  stable  range  of  emitted  currents. 

In  Fig.  (4-2a)  the  emitted  current  from  a  stable  meniscus  is  shown. 

The  current  waveform  was  obtained  by  monitoring  the  particles  impinging  on 

the  T0F  collector  by  means  of  an  oscilloscope.  Photomicrographs  of  the 

capillary  and  stable  meniscus  arc-  shown  in  Fig.  (4-2b)  .  Apex  angles 

ranging  from  60  to  110  degrees  were  detected.  At  low  current  levels 

(approximately  50  microamps)  conical  meniscuses  with  apex  angles  very  near 

o 

to  Taylor's  predicted  value  of  98.6  were  observed.  At  higher  current 
levels  the  meniscus  usually  had  a  smaller  apex  8ngle. 
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4 . 2  Capillary  Current  Measurement 

The  capillary  current  was  measured  with  a  0-500  microamp  meter  placed 
near  the  capillary  to  reduce  the  possibility  of  leakage  currents.  The 
field  evaporation  theory  developed  in  Chapter  2  predicts  that  the  log¬ 
arithm  of  the  emitted  current  varies  as  the  square  root  of  the  applied 
field.  To  verify  this  relationship  it  was  necessary  to  obtain  the  func¬ 
tional  dependence  of  the  emitted  current  on  the  applied  voltage.  The 
current-voltage  (I-V)  characteristic  was  experimentally  determined  by 
measuring  the  current  and  voltage  at  10  microamp  intervals  and  is  plotted 
in  Fig.  (4-3)  for  currents  ranging  from  20  to  500  microamps.  By  plotting 
the  logarithm  of  the  current  as  a  function  of  the  square  root  of  the 
voltage  the  predicted  theoretical  relationship  is  a  straight  line. 

The  experimental  I-V  relationship  is  divided  into  three  regions: 

(1)  Space  charge  limited:  (2)  Emission  Limited:  and  (3)  Saturation  limited. 
At  low  currents  and  low  accelerating  voltages  the  ion  emission  is  space 
charge  limited  and  the  I-V  characteristic  approximately  follows  the  three- 
halves  power  law.  At  100  microamps,  the  capillary  current  enters  an 
emission  limited  region.  In  this  mode  of  operation  the  I-V  characteristic 
given  by  a  1 ea st -squa res-error  fit  to  the  data  points  is 


-30  1/2 

I  =  0.99  x  10  EXPC0.75  V  ) 


(4-1) 


In  the  range  of  100  to  300  microamps  the  number  of  lattice  sites 
undergoing  evaporation  remains  essentially  constant  and  the  increase  in 
ion  current  is  produced  by  decreasing  the  time  constant  for  evaporation 
at  those  sites.  Hence,  the  nearly  exponential  I-V  relationship  results 
Above  300  microamps  the  source  begins  to  saturate.  Additional 
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increases  in  beam  current  require  additional  sites  to  undergo  evaporation. 
As  the  number  of  evaporating  sites  increases  the  radius  of  curvature  at 
the  new  sites  increases',  and  hence,  the  voltage  necessary  to  produce 
the  required  evaporating  field  strength  increases.  This  results  in  a 
dcviatior  in  the  I-V  curve  from  the  exponential  at  high  current  levels. 

By  considering  the  exponential  region  of  the  I-V  curve,  the  radius 
and  binding  energy  of  the  emitting  surface  can  be  obtained.  After  in¬ 
serting  the  numerical  values  for  the  physical  constants  at  300°  K,  the 
theoretical  I-V  dependence  given  by  Eqn.  (2-24)  is 

I  =  1.6  x  l(f6  N  EXP(-40Qo)  EXP(1.5  x  10_3£1/2>  (4-2) 

By  equating  the  coefficients  of  the  exponentials  of  Eqns.  (4-1)  and 
(4-2),  the  binding  energy  at  the  liquid  surface  can  be  obtained. 

_  OQ  _  g 

0.99  x  10  =  (1.6  x  10  )  N  EXP( -40Q  )  (4-3) 


In  order  to  solve  for  the  binding  energy,  Q^,  the  number  of  lattice 

sites  undergoing  evaporation,  N,  must  be  determined.  For  an  exponential 

I-V  characteristic,  N  must  be  constant  throughout  the  range  of  operation 

The  maximum  current  that  can  be  emitted  from  a  single  latticesite  is 

-13 

given  by  Eqn.  (2-23)  assuming  T  =10  seconds. 


=  1.6  x  10 


-6 

amp 


(4-4) 


At  300  ft  A  (the  maximum  current  in  the  exponential  region  of  operation) 
the  minimum  number  of  lattice  sites  is  188.  Using  this  value  for  N 
Eqn.  (4-3)  can  be  solved  for  the  binding  energy. 


Q 


o 


1.52  ev  . 


(4-5) 
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An  approximate  value  lor  the  radius  at  the  tip  of  the  liquid  cone  can 
be  obtained  by  equating  the  coefficients  of  the  arguments  of  the  expo¬ 
nentials  in  Eqns.  (4-1)  and  (4-2). 

1.5  x  10‘3£  1/2  =  0.75  V1/2  (4-6) 

In  order  to  solve  this  equation,  the  relationship  between  the  applied 
voltage  and  electric  field  at  the  tip  of  the  liquid  cone  given  by  Eqn. 

(2-13)  must  be  used.  Eqn.  (2-13)  was  derived  for  a  system  in  which  the 
cathode  was  an  equipotent lal  surface  in  the  conical  geometry.  In  the 
experimental  geometry,  the  planer  accelerating  electrode  does  not  possess 
the  shape  of  an  equipotential  in  conical  geometry;  therefore,  the  inter¬ 
electrode  spacing  a^,  cannot  be  directly  determined.  Intuitively,  one 
would  surmise  that  the  effect  of  the  planer  cathode  and  the  cylindrical 
capillary  would  be  to  concentrate  the  field  lines  at  the  conical  tip  and 
thus  enhance  the  field. 

An  approximate  value  for  the  electrode  spacing  can  be  obtained  by 

determining  the  minimum  applied  voltage  that  is  necessary  to  support  the 

conical  geometry.  For  the  gal 1 ium- indium  alloy  the  voltage  is  given 

(27) 

in  terms  of  the  electrode  spacing  by 

V  -  3.58  x  1 04  (a,)12  (4-7) 

o  1 

The  minimum  voltage  needed  for  the  conical  geometry  was  experimentally 
found  to  be  6.35  KV.  This  corresponds  to  an  eflective  i nterel ect rode 
spacing  of  approximately  =  0.031  cm.  whereas  the  actual  geometric 
separation  was  0.15  cm.  Therefore,  the  distorted  electrode  geometry 


enhanced  the  field  ai  the  tip  of  the  cone. 


61 


Using  the  ab/'c  value  for  a  ,  Eqn .  (2-13)  can  be  used  to  determine 
the  electric  field-voltage  relationship  at  the  tip  ol  the  cone.  Further, 
by  substituting  Eqn.  (2-13)  into  (4-6)  the  radius  at  the  tip  of  the  liquid 
cone  was  determined  to  be 


Po 


=  520  X 


(i-8) 


The  upper  and  lowe*"  limits  on  the  range  of  the  emitted  current  were 
dictated  by  the  corresponding  ranges  of  the  applied  accelerating  voltage. 

At  capillary  currents  in  the  20  to  50  microamp  range  the  applied  voltage 
was  marginal  with  respect  to  support  of  the  conical  geometry.  At  these 
low  current  levels  the  meniscus  would  at  times  oscillate  due  to  a  Hayleigh- 
Taylor  instability  as  given  by  Eqn.  (1-21).  The  lowest  unstable  mode  for 
a  liquid  hemisphere  at  a  capillary  tip  is  the  n  =  2  mode.  For  the  Gallium- 
indium  working  fluid  and  a  125  f»m  radius  capillary  the  unstable  voltage 
is  given  by 

1  /? 

V  =  (5  'Y  r.a  (300 )'2 )  =  6.5  KV  (4-9) 

This  value  of  voltage  for  the  n  =  2  instability  is  within  3%  of  the 
6.31  to  6.35  KV  range  of  applied  voltages  needed  to  produce  20  -  50 
microamps  of  emitted  current.  (See  Fig.  (4-3)). 

In  this  mode  of  operation  the  emitted  current  was  observed  as  a 
train  of  periodic  pulses  as  shown  in  Fig.  (4-4a)  for  a  current  of  20j<A. 

The  expansion  of  a  single  pulse  in  Fig.  (4-4b)  shows  fiat  each  pulse  was 
ac'ually  composed  of  mBny  charged  particles.  As  the  current  level  was 
changed  the  pulses  remained  nearly  constant,  but  the  frequency  of  the 
pulses  varied.  Photomicrographs  of  the  liquid  meniscus  shown  in  Fig.  (4-4c) 
depict  the  oscillatory  motion  of  the  n  =  2  instability. 


FIG.  4-4  LOW  VOLTAGE  INSTABILITY,  (a)  EMITTED  CURRENT 
PULSES  AT  1  »s/cm  SWEEP  SPEED,  (b)  SINGLE 
PULSE  AT  20  fis/cm  SWEEP  SPEED  'c)  PHOTO¬ 
MICROGRAPHS  OF  MENISCUS  AT  2  '  INTERVALS. 


At  currents  in  excess  of  500  microamps  the  second  unstable  Ruvleigh- 
Ta  ylor  mode  produced  oscillations  in  the  liquid  cone.  In  the  n  =  2  mode, 
the  liquid  cone  underwent  longitudinal  oscillations.  Fur  the  second 
unstable  mode,  n  —  3,  the  oscillations  were  in  a  transverse  direction. 

In  this  mode  the  unstable  applied  voltage  is  given  by 

4  1/2 

V  9  x  10  )  =  7.7  KV  (4-10) 

Experimentally  this  instability  was  detected  at  approximately 
6.88  KV.  The  predicted  unstable  voltage  is  higher  than  the  experimental 
value  by  10%.  This  discrepancy  cun  be  explained  by  noting  that  in  the 
derivation  of  Eqn.  (1-21)  the  electric  field  normal  to  the  meniscus  w is 
assumed  to  be  equal  to  the  field  about  a  sphere.  The  effect  of  the 
conical  shape  would  be  to  concentrate  the  field  lines  st  the  tip  of  the 
hemisphere  and  thus  lower  the  unstable  voltage.  Photomicrographs  oi  the 
n  =  3  instability  at  a  capillary  current  oi  750  microamps  are  shown  in 
Fig.  (4-5). 

4 . 3  Specific  Charge  Measurement 

The  specific  charge  distribution  of  the  field  evaporated  particles 
was  determined  by  use  of  both  the  t ime-of -f 1 lght  mass  specti ometer  (T0FMS) 
and  the  quadrupole  mass  filter.  Data  points  were  obtained  at  50  microsmp 
intervals  for  capillary  currents  ranging  from  50  to  500  microamps.  Over 
this  entire  range  of  operating  currents  the  specific  charge  distribution 
remained  constant  to  within  the  accuracy  of  the  detection  apparatus,  which 
was  approximately  l%of  the  total  detected  current. 

The  quadrupole  mess  filter  was  used  to  obtain  high  resolution  spectra 
of  near  axial  particles  and  thus  sampled  only  a  very  smjll  portion  of  the 
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beam.  A  typical  specific  charge  distribution  is  shown  in  Fig.  (4-6).  The 
spectrum  was  found  to  consist  of  75%  gallium  dimers,  24%  indium  dimers  and 
less  than  1%  gallium  tetramers.  The  corresponding  specific  charges  are 
shown  in  Table  (4-1). 


TABLE  (4-1) 

SPECIFIC  CHARGE 

;  DISTRIBUTION 

Ionic  Species 

Percentage 

Specific  Charge 

»,* 

75 

7.1  x  10s  coul/kg 

In2  + 

24 

4.3  x  10s  coul/kg 

1 

3.5  x  10s  coul/kg 

The  specific  charges  of  the  field  evaporated  particles  were  all 
found  to  lie  between  (3.5  -  7.1)  x  105  coul/kg.  The  average  specific 
charge  of  the  beam  was  6.3  x  10s  coul/kg. 

The  quadrupole  mass  filter  data  was  corroborated  by  the  TOFMS.  The 
TOF  device  was  capable  of  detecting  all  particles  spraying  within  a  cone 
of  30°  apex  angle,  but  it  was  inherently  a  low  resolution  device.  A  typical 
TOF  collector  waveform  is  shown  in  Fig.  (4-7a).  By  numerically  differen¬ 
tiating  this  curve,  the  specific  charge  distribution  is  /otained  as  per 
Eqn .  (3-11)  and  is  plotted  in  Fig.  (<l-7b)  . 

The  TOF  spectrum  has  broad  peaks  centered  at  5.1  x  1CS  coul/kg  and 
8.1  x  10s  coul/kg,  with  amplitudes  of  22%  and  78%  of  the  total  current 
respectively.  The  width  of  these  peaks  and  their  deviations  from  the 
quadrupole  data  can  be  explained  by  considering  the  basic  limitations  of 
the  TOF  measurement  technique.  First  looking  at  Fig.  (4-7a),  the  collector 
waveform  is  not  rectangular  (as  it  /ould  be  ideally),  but  it  is  rounded. 

This  rounding  tends  to  obscure  the  peaks  of  the  spectrum  and  broaden  them. 


It  was  caused  b.  the  finite  rise  time  of  the  collector  circuitry,  and  the 
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FIG.  4-6  SPECIFIC  CHARGE  DISTRIBUTION  DETERMINED 
BY  QUADRUPOLE  MASS  FILTER 
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high  velocity  of  the  emitted  particles  which  necessitated  a  0.5  micro¬ 
second/cm  sweep  speed.  The  resolution  was  further  limited  by  the  line 
width  of  the  oscilloscope  trace. 

The  second  fundamental  limitation  is  the  establishment  of  a  time 
base  for  the  oscilloscope  sweep  relative  to  the  high  voltage  pulse.  It 
is  necessary  for  the  high  voltage  to  start  to  decline  in  order  to  trigger 
the  oscilloscope  sweep.  Thus  the  zero  on  the  oscilloscope  display  corres¬ 
ponds  to  some  positive  instant  in  time  relative  to  the  high  voltage  pulse, 
and  the  entire  waveform  fs  shifted  in  time.  A  shift  in  the  time  scale  of 
200  nanoseconds  would  align  the  TOF  spectra  with  the  higher  resolution 
quadrupole  data. 

The  limitation  in  resolution  is  justified  by  the  large  cross-section 
of  the  beam  that  the  TOFMS  can  detect.  Since  the  composition  of  the 
working  fluid  is  known,  only  gallium  ana  indium  ions  can  be  present  in 
the  beam.  Thus,  the  accuracy  of  the  TOF  technique  is  sufficient  to  deter¬ 
mine  the  particular  specilic  charges  in  the  beam,  and  it  can  be  assumed 
that  the  two  TOF  peaks  correspond  to  gallium  and  indium  dimers. 


4 .4  Formation  Energy  Measurement 

The  formation  energy  of  the  field  evaporated  ions  was  determined  by 
utilizing  the  retarding  potential  cell.  Data  points  were  taken  at  10 
volt  intervals  in  the  retarding  voltage  while  maintaining  the  accelerating 
voltage  constant  to  within  10  volts.  Data  runs  were  made  at  50,  100,  200, 
300,  400,  and  500  microamps  capillary  current.  At  all  current  levels  the 
formation  energy  distribution  was  approximately  gaussian.  As  the  current 
was  increased  from  50  to  200  microamps,  the  position  of  the  maximum  and 


the  variance  of  the  distribution  increased.  In  the  200  to  500  microamp 
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range  the  energy  distribution  remained  constant  to  within  the  10  volt 
accuracy  of  the  measurement  technique. 

The  retarding  potential  current-voltage  characteristics  are  shown  in 
Fig.  (4-8).  Dy  differentiating  these  curves  the  energy  distributions 
are  obtained  and  they  are  plotted  in  Fig.  (4-9).  The  maximum  of  the 
energy  distributions  and  the  corresponding  variances  are  listed  in 
Table  (4-2)  . 

TABLE  (4-2)  FORMATION  ENERGY  DISTRIBUTION 

Current  Peak  Energy  Variance 

2 

50  pA  60  volts  (38) 

o 

100  75  volts  (42)* 

(47)2 


200-500  uA 


85  volts 
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£  LERATING  POTENTIAL  MINUS  RETARDING  POTENTIAL  IN  VOLTS 

FIG.  4-8  RETARDING  POTENTIAL  CURRENT-VOLTAGE  CHARACTERISTIC 


FIG.  4-9  ENERGY  DISTRIBUTION  OF  FIELD  EVAPORATED  PARTICLES 


5.  CONCLUSIONS 
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This  study  has  shown  that  under  the  appropriate  operating  conditions, 
EHD  spraying  techniques  can  be  utilized  for  the  production  of  high 
current  density  ion  beams.  The  effects  of  capillary  wetting  and  the 
feedback  control  (current  regulation)  of  the  accelerating  voltage  were 
responsible  for  the  performance  obtained.  Without  these  techniques  ions 
and  heavy  particles  would  have  been  simultaneously  produced,  and  the 
present  wide  dynamic  range  of  operation  currents  would  not  have  been 
attained . 

5.1  Comparison  with  Previous  Results 

Several  close  analogies  can  be  drawn  between  the  operation  of  an 

EHD  spray  operating  in  the  field  emission  mode  and  the  results  of  previous 

experiments  utilizing  an  EHD  spray  for  the  production  of  heavy  charged 

particles.  The  instabilities  of  the  liquid  meniscus  during  both  low  and 

high  field  operation  correlate  very  closely  with  the  extension  of  Rayleigh's 

(18) 

theory  derived  by  Hendricks  et  si.  The  close  experimental  agreement 

(3%)  with  theory  for  T’h°  minimum  spraying  potential  for  the  n  =  2 
instability  complements  previous  results  using  liquids  with  low  conducti¬ 
vity.  Hence,  for  a  range  of  15  orders  of  magnitude  of  conductivities, 
the  modification  of  Rayleigh's  theory  applicable  to  a  liquid  hemisphere 
at  the  tip  of  an  electrified  capillary  has  been  experimentally  verified. 

Thus  except  for  the  number  of  unstable  eigen-frequencies  allowed,  the 
behavior  of  an  electrodynamically  unstable  liquid  hemisphere  appears  to 
be  physically  the  same  as  that  of  an  electrodynamically  unstable  free 
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charged  sphere. 

The  operution  of  an  EHD  spray  at  voltages  near  the  minimum  spraying 

( 6  ) 

potential  was  previously  observed  by  Carson.  '  The  longitudinal 

oscillation  of  the  liquid  meniscus  at  low  voltages  corresponds  to  Carson's 

"naturally  pulsed"  mode  of  operation.  Oscilloscope  traces  of  the  emitted 

current  pulses  have  indicated  that  each  pulse  is  actually  composed  of 

many  charged  particles.  This  result  is  also  in  agreement  with  Carson  but 

(37) 

contradicts  Pierce's  theory  that  only  one  particle  is  produced  per 

pulse-  This  may  be  the  result  of  an  inability  to  resolve  the  major  current 

pulses  into  a  fine  structure  which  would  support  Pierce's  contention. 

The  physical  phenomenon  responsible  for  the  production  of  these 

pulsations  has  not  been  determined.  Several  possible  mechanisms  in  low 

conductivity  liquids  fi.e.  charge  relaxation  time  and  thermal  effects) 

(38) 

have  been  studied  by  Hendricks  and  Crowley.  The  observed  pulsations 
with  the  gallium-indium  working  fluid  precludes  the  effects  of  electrical 
and  thermal  time  constants  in  the  liquid.  It  may  be  that  the  frequency 
of  the  observed  pulses  for  high  conductivity  liquids  is  governed  by  the 
external  time  constant  of  the  accelerating  voltage  supply. 

Unlike  the  EHD  spraying  of  low  conductivity  liquids  where  the  charged 
particle  beam  has  been  assumed  to  be  produced  from  an  unstable  meniscus, 
the  emission  of  ions  from  the  gallium-indium  working  fluid  was  by  a 
direct  interaction  of  the  applied  field  on  the  surface  without  a  bulk 
mechanical  oscillation.  Hence,  a  different  functional  relationship 
between  the  specific  charge  of  the  emitted  particles  and  the  applied 

(19) 

field  was  obtained.  Pfiefer's  energy  minimization  analysis  of  an 

unstable  liquid  plane  predicts  a  specific  charge  varying  as  the  one-hslf 
power  of  the  applied  field  as  shown  in  Eqn.  (1-22).  As  the  field  intensity 
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is  increased  the  current  is  increased  due  to  an  increase  in  the  specific 
charge  of  the  emitted  particles.  In  the  field  evaporation  mode  of 
spraying  the  specific  charge  of  the  emitted  particles  is  constant,  and 
the  current  is  increased  by  increasing  the  number  of  particles  produced 
per  second. 

5 . 2  Critique  of  Theoretical  Considerations 

In  this  section  the  models  developed  to  explain  the  field  evaporation 

mode  of  EHD  spraying  will  be  summarized.  The  predicted  conical  meniscus 

was  observed  for  a  wide  range  of  applied  electric  field  intensities. 

Taylor's  theoretical  apex  angle  of  98.6°  was  observed  along  with  angles 
o  o 

ranging  from  60  to  110  .  In  general  the  angle  appeared  to  be  determined 
by  the  mass  of  working  fluid  in  the  meniscus  and  not  by  the  Taylor 
stability  criteria.  The  present  experiment  differed  from  that  assumed 
in  Taylor's  derivation  in  three  important  aspects:  (1)  The  field  about 
the  cone  was  not  uniform;  (2)  The  meniscus  was  not  a  free  surface  but  was 
constrained  by  the  capillary;  and  (3)  The  internal  pressure  within  the 
meniscus  was  non-zero.  These  three  effects  combined  to  allow  the  exper¬ 
imentally  observed  range  of  "stable"  apex  angles. 

The  experimental  current-voltage  relationship  given  in  Eqn .  (4-1) 
agrees  with  the  functional  dependence  predicted  by  the  field  evaporation 
theory.  The  exponential  nature  of  the  I-V  curve  tends  to  make  curve 
fitting  very  inaccurate;  therefore,  the  experimental  I-V  curve  in  itself 
does  not  justify  the  use  of  the  field  evaporation  model.  Further  evidence, 
however,  supporting  the  field  evaporation  theory  is  gained  from  the  specific 
charge  distribution  of  the  emitted  particles  and  the  electric  field  inten¬ 
sity  at  the  conical  surface.  The  ion-sized  particles  produced,  the 
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extremely  high  field  intensities  measured,  and  the  experimental  I-V  curve 
all  serve  to  indicate  that  a  field  evaporation  phenomenon  is  responsible 
for  the  production  of  the  high  current  density  ion  beam. 

The  electric  field  intensity  at  the  apex  of  the  meniscus  can  be 
calculated  from  Eqn.  (2—13)  by  assuming  the  radius  at  the  apex  of  the  cone 
to  be  520  X  as  determined  in  Eqn.  (4-8). 

£  =  2.5  x  10*  V  (5-1) 

0 

The  electric  field  intensity  for  a  current  of  300  is  then  1.67  x  10 
volt  s/meter . 

Over  forty  elements  have  been  field  evaporated^30 ^  In  Fig.  (5-1), 
the  evaporating  field  strengths  of  a  representative  group  of  these 
elements  are  plotted  as  a  function  of  their  binding  energies.  The  present 
data  for  the  gallium-indium  working  fluid  is  included  using  the  experi¬ 
mentally  determined  values  for  the  electric  field  intensity  and  binding 
i  energy. 

The  only  significant  loss  mechanism  in  the  ion  beam  formation  was 
through  the  space  charge  of  the  emitted  particles.  Losses  due  to  both 
the  binding  energy  of  the  atoms  and  the  finite  conductivity  of  the  working 
fluid  were  found  to  be  negligible  in  comparison  to  the  kinetic  energy  of 
the  ions.  The  predicted  potential  loss  due  to  the  space  charge,  as 
given  by  Eqn.  (2-38)  and  the  peak  values  of  the  experimental  retarding 
potential  curves  are  listed  in  Table  (5-1)  as  a  function  of  the  emitted 


current . 
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TABLE  (5-1) 

THEORETICAL  AND  EXPERIMENTAL  FORMATION  POTENTIALS 

Theoret.  ical 

Experimental 

I  CiiA) 

Fotantial  (V) 

Potential  (V) 

50 

16 

60 

100 

31 

75 

200 

62 

85 

300 

9-> 

85 

400 

121 

85 

500 

151 

85 

The  poor  agreement  between  the  theoretical  and  experimental  data 
indicates  that  the  proposed  space  charge  model  is  at  best  a  poor  appro¬ 
ximation  to  the  actual  physical  phenomenon.  The  basic  assumptions  in 
the  derivation  of  Eqn.  (2-38)  ol  a  constant  current  density  and  a  parallel 
ion  beam  were  violated  by  the  experimental  system  where,  in  essence,  the 
beam  was  diverging  from  a  near  point  source.  Also,  the  total  emitted 
current  and  the  current  density  were  not  linearly  related.  As  the  current 
was  increased  the  emitting  area  also  increased. 


5 . 3  Thruster  Utrliza r ion 

The  usefulness  of  the  field  evaporated  dimer  beam  as  a  microthruster 
can  be  determined  by  considering  the  following  parameters:  (1)  Thrust; 
(2)  Parti  cl e  formation  efficiency;  and  (3)  Beam  conversion  efficiency. 

An  electrostatically  accelerated  beam  of  gallium  and  indium  dimers 
produces  a  thrust  given  by  Eqn.  (1-4). 


-3 


IV 


a 


F  =  1.78  x  10 


1/2 


(5-2 
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The  20  -  500  microamp  range  of  emitted  currents  that  can  be  produced  leads 
to  a  range  of  thrusts  varying  from  2.8  to  73.8  micronewtons  (1.3  to  33.5 
micropounds).  The  thrust  is  listed  as  a  function  of  the  current  and 
accelerating  voltage  in  Table  (5-2). 

TABLE  (5-2)  THRUST  AS  A  FUNCTI0N~0F  CURRENT  AND  ACCELERATING  VOLTAGE 


(MA) 

V  (kv) 

F  (^N) 

20 

6.31 

2.8 

50 

6.35 

7.1 

100 

6.44 

14.4 

200 

6.57 

28.9 

300 

6.68 

43.6 

400 

6.78 

58.6 

500 

6.88 

73.8 

As  the  thrust  is  increased  the  accelerating  voltage  remains  nearly 
constant;  therefore,  the  velocity  of  the  emitted  particles  remains 
essentially  constant,  and  they  could  be  maintained  at  a  near  optimum 
velocity  independent  of  the  thrust  level. 

The  particle  formation  efficiency  of  the  charged  particle  beam  can 
by  computed  by  substituting  the  retarding  potential  data  into  Eqn.  (1-5) 
The  formation  efficiency  is  listed  as  a  function  of  the  emitted  current 


in  Table  (5-3) . 
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TABLE  (5-3)  FORMATION  EFFICIENCY  AS  A 

FUNCTION  OF  EMITTED  CURRENT 

l  Cm A) 

T){  (%) 

50 

99.1 

100 

98.8 

200 

98.7 

300 

98.8 

400 

98.8 

500 

98.9 

The  approximately  1%  loss  in  efficiency  was  previously  attributed  to 
space  charge  effects.  This  energy  loss  appears  as  a  transverse  accelera¬ 
tion  of  the  ions  and  leads  to  beam  spreading.  Such  a  divergent  beam  has 
reduced  thrust  producing  capabilities  compared  to  a  parallel  beam.  It 
should  also  be  noted  here  that  all  of  the  ions  were  not  emitted  with 
parallel  trajectories.  In  addition  to  the  losses  given  above  there  was 
Bn  inherent  beam  spreading  due  to  the  electric  field  geometry  at  the 
capillary  tip. 

The  constant  specific  charge  distribution  of  the  particles  produced 
by  the  field  evaporation  process  leads  to  a  constant  beam  conversion 
efficiency  independent  of  the  thrust  level.  By  substituting  the  experi¬ 
mentally  determined  specific  charge  distribution  into  Eqn.  (1-12),  the 
beam  conversion  efficiency  was  found  to  be 

=  97.5%  (5-3) 

b 

Thus,  the  distribution  of  specific  charges  present  in  the  beam  degrades 
the  thrust  producing  capability  of  the  system  by  2.5%. 

In  summary,  the  ion  beam  produced  by  EHD  spraying  techniques  has 
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been  shown  to  be  useful  as  a  source  of  ionic  particles  for  micro-thruster 
utilization.  The  high  thrust  produced,  high  formation  efficiency,  and 
high  beam  conversion  efficiency  lead  to  a  system  combining  the  advantages 
of  existing  ionic  and  heavy  particle  electrostatic  thrusters. 

5 . 4  Recommendations  for  Future  Research 

The  present  solutions  to  the  problems  encountered  in  the  course  of 
this  study  have  left  many  questions  unanswered  and  have  produced  many 
more  new  areas  of  interest  for  future  research.  Some  of  the  more  pertinent 
problems  are  described  below. 

(1)  Apply  feedbcek  r  ^ntrol  (current  regulation)  to  the  EHD  spraying 
of  various  working  fluids.  Construct  conical  electrodes  so  that  a  uniform 
field  is  obtained  and  TayJor's  stability  theory  can  be  verified.  Analyze 
the  parameters  of  the  sprayed  particles  as  a  function  of  liquid  conducti¬ 
vity,  surface  tension  and  viscosity. 

(2)  Obtain  a  solution  of  Child's  law  for  the  conical  electrode 
configuration  so  that  a  precise  analysis  of  space  charge  effects  may  be 
carried  out . 

(3)  Extend  Taylor's  conical  stability  criterion  from  the  free  surface 
case  to  the  case  of  a  liquid  meniscus  constrained  by  a  capillary. 

(4)  Using  the  present  experimental  geometry,  study  the  field  evapora¬ 
tion  of  other  conducting  liquids  such  as  mercury,  cesium,  sodium,  and 
potassium  to  determine  their  usefulness  in  electrostatic  thrusters. 

(5)  Spray  pure  gallium  instead  of  the  alloy  and  determine  the  effect 
of  temperature  on  the  field  evaporation  process. 

(6)  Explore  the  possibility  of  producing  higher  currents  and  thrusts 
by  utilizing  a  different  source  geometry  such  as  an  annulus,  thin  slit, 
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or  rotating  disc. 

(7)  The  addition  of  an  AC  component  at  an  unstable  Rayleigh  eigen- 
frequency  to  the  accelerating  voltage  could  enhance  the  breakup  of  the 
meniscus  into  larger  particles.  Such  particles  could  be  useful  in  heavy 
particle  propulsion  systems. 

(8)  Study  post  acceleration  or  deceleration  of  the  field  evaporated 
ions  as  a  means  of  obtaining  an  optimum  exhaust  velocity. 

(9)  Determine  the  effects  of  intense  negative  fields  on  liquid 
metals.  Preliminary  studies  produced  a  10  mA  electron  source  from  a 
single  capillary. 

(10)  Analyze  the  effect  of  the  accelerating  voltage  supply  impedance 
on  the  "naturally  pulsed"  mode  of  EHD  spraying  for  both  conducting  and 
non-conduct  ng  liquids. 

(lil  Study  the  systems  aspects  of  microthrusters.  Construct  a 
■working  thruster  module  composed  of  an  array  of  capillaries  operating  in 
the  field  evaporation  mode  capable  of  producing  0.5  millipounds  of  thrust 
Perform  10,000  hour  life  tests  to  determine  reliability  for  space  applies 
tions.  Analyze  the  compatibility  of  above  system  with  respect  to  thrust 
vectoring,  propellant  feed  systems,  and  power  conditioning  of  the  output 
from  present  photovoltat ic  sources. 


(39) 
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